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Abstract 

Delaunay flip is an elegant, simple tool to convert a triangulation of a point set to its 
Delaunay triangulation. The technique has been researched extensively for full dimensional 
triangulations of point sets. However, an important case of triangulations which are not 
full dimensional is surface triangulations in three dimensions. In this paper we address the 
question of converting a surface triangulation to a subcomplex of the Delaunay triangulation 
with edge flips. We show that the surface triangulations which closely approximate a smooth 
surface with uniform density can be transformed to a Delaunay triangulation with a simple 
edge flip algorithm. The condition on uniformity becomes less stringent with increasing 
density of the triangulation. If the condition is dropped completely, the flip algorithm still 
terminates although the output surface triangulation becomes "almost Delaunay" instead 
of exactly Delaunay. 
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1 Introduction 



The importance of computing Delaunay triangulations of point sets in applications of science 
and engineering cannot be overemphasized. A number of different algorithms have been pro- 
posed for their computations O [HI [H] . Among them flip based algorithms are most popular 
and perhaps the most dominant approach in practice. The sheer elegance and simplicity of this 
approach make it attractive to implement. 

Since the introduction of Delaunay flips by Lawson |17j for point sets in two dimensions, 
several important extensions have been made including higher dimensional point sets. Given 
any triangulation of the convex hull of a point set in two dimensions, it is known that Delaunay 
edge flips can convert the triangulation to the Delaunay triangulation. The rule for Delaunay 
edge flips is simple. First, check locally if the circumscribing ball of a triangle t contains a 
vertex of another triangle t' sharing an edge e with it. If so, replace e with the other diagonal 
edge contained in the union of t and t' . An elegant result is that this process terminates with 
the output as the Delaunay triangulation [5l [HI [ll]. In higher dimensions, the edge flips can 
be naturally extended to bi-stellar flips [13] • However, the approach extended by bi-stellar 
flips does not work in higher dimensions. Already in three dimensions there are examples 
where the flipping process can get stuck [T3]. Notwithstanding this negative result, Joe [16] 
showed how to construct Delaunay triangulations by inserting points incrementally and applying 
bi-stellar flips after each point insertion. Edelsbrunner and Shah [15] extended this result to 
higher dimensional point sets and to weighted Delaunay triangulations. Recently, Shewchuk |19j 
showed that a combination of flips and some other local operation called star splay can convert 
an "almost Delaunay" triangulation to the Delaunay one quite efficiently. 

All the aforementioned results deal with full dimensional triangulations of a point set. An 
important case of triangulations which are not full dimensional is surface triangulations in 
three dimensions. Given the increasing demand of computing surface triangulations that are 
sub-complexes of Delaunay triangulations [H [71 [12] , it is natural to ask if a surface triangulation 
can be converted to a Delaunay one by edge flips and, if so, under what conditions. Such a flip 
algorithm will be useful in many applications. For example, in geometric modeling, shapes are 
often represented with subdivision surfaces [20] or with isosurfaces jl8j . These surfaces are not 
necessarily Delaunay. If one can convert these surfaces to a Delaunay one, a number of tools 
that exploit Delaunay properties can be used for further processing. 

In this work we address the question of Delaunay flips in surface triangulations. Notice that 
our goal is to convert a surface triangulation embedded in to another surface triangulation 
which is a sub-complex of the three dimensional Delaunay triangulation of the vertex set. This is 
different from the framework considered elsewhere [6l [13] where a triangulated surface endowed 
with a flat metric is converted into an intrinsic Delaunay triangulation comprised of simplices 
(not necessarily planar) embedded in the surface. In this case the embedding of the surface in 

does not play any role whereas in our case the positions of the vertices in R'^ determine the 
Delaunay flips. 

It turns out that dense surface triangulations are amenable to a simple edge flip algorithm. 
A triangulation is dense if it approximates a smooth surface with sufficient resolution. We show 
that a dense triangulation can be ffipped to a Delaunay triangulation if the density is uniform 
in some sense. The condition on uniformity depends on the density. The higher the density, the 
less stringent is the condition. The practical implication of this result is that reasonably dense 
triangulations can be converted to Delaunay triangulations with a simple edge flip algorithm. 
Such dense triangulations are numerous in practice. Subdivision and isosurface meshes are two 
such examples. Triangulations of moving vertices offer another such example [10\ I19j. In fact. 
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the results in this paper have been used for a recent algorithm on maintaining deforming meshes 
with provable guarantees [10] . What happens if we do not have the uniformity condition? We 
show that the flip algorithm still terminates but the output surface may not be Delaunay. 
Nonetheless, this surface is "almost Delaunay" in the sense that the diametric ball of each 
triangle shrunk by a small amount remains empty. Bandyopadhyay and Snoeyink [3] showed the 
usefulness of such approximate Delaunay triangulations in molecular modeling. Because of the 
approximate emptiness properties of the circumscribing balls of the triangles, these approximate 
Delaunay triangulations may find other applications where exact Delaunay triangulations are 
not required. 

2 Preliminaries 

2.1 Definitions and results 

We need some definitions and results from e-sampling theory [1]. Let S C be a smooth 
compact surface without boundary. The medial axis is the set of centers of all maximally 
empty balls. The reach 7(S) of S is the infimum over Euclidean distances of all points in T, to 
its medial axis. This is also the infimum of the local feature size of S as defined by Amenta 
and Bern [1]. 

A surface triangulation T is a finite simplicial 2- manifold embedded in M^. We say T is a 
triangulation of a surface S if vertices of T lie in S and its underlying space |T| is homeomorphic 
to S. The triangulation T has a consistent orientation with S if the oriented normal of each 
triangle makes at most ^ angle with the oriented normals of S at the vertices. For a triangle 
t £ T, let p{t) denote the circumradius of t. 

2.2 Uniform dense triangulations 

Definition 1 A triangulation T of a surface E is e-dense if each triangle t £ T has p(t) ^ e7(S) 
and T has a consistent orientation with S. Furthermore, for 6 < 1, if any two vertices in T 
has distance more than 5e^{T,), T is called {e, 5)-dense or 6-uniform and e-dense. 

We use notation B(c,r) to denote a ball with center c and radius r. A circumscribing ball 
of a triangle t C is any ball that has the vertices of t on its boundary. The diametric ball 
Dt is the smallest such ball; Dt = B{c,p(t)) where c is the circumcenter of t. We say a vertex 
V (zT stabs a ball B if v lies inside B. If t shares an edge, say pq, with a triangle t' = pqs, then 
s is a neighbor vertex of t. Clearly, each triangle has three neighbor vertices. 

Definition 2 A triangle t £ T is stabbed if Dt is stabbed by a vertex ofT. We say t is locally 
stabbed if the stabbing vertex is one of the three neighbor vertices oft (FigureU^. 

Theorem 12.11 and Theorem 12.21 are main results for uniform dense triangulations. 

Theorem 2.1 For 5 = 2sin24e and e < any (e, 6)-dense triangulation has a stabbed trian- 
gle if and only if it has a locally stabbed triangle. 

Notice that the bound on 6 is 0(e). This implies that dense triangulations require only mild 
restrictions on its uniformity. 

We will flip edges of dense triangulations to make it Delaunay. Suppose pq is an edge in a 
surface triangulation T. Flipping pq means replacing two triangles, say pgr and pqs, incident 
to pq in T by the triangles prs and qrs. If the new triangulation is T' we write T ^ T' . 
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A triangulation T is flipped to a triangulation T' if there is a sequence of edge flips so that 

Definition 3 An edge in a surface triangulation is called flippable if it is incident to a locally 
stabbed triangle. 

One can devise an easy algorithm to convert a (e, (5)-dense triangulation to a Delaunay 
triangulation using Theorem l2.11 Simply flip any flippable edge existing in the triangulation. If 
one can prove that this flip algorithm terminates and (e, 5)-density is maintained after each flip, 
we will have an algorithm to flip an (e, (5)-dense triangulation to a triangulation that does not 
have any stabbed triangle. This means each triangle in the new triangulation has its diametric 
ball that does not contain any vertices of T. In other words, this new triangulation is Delaunay. 
Actually, the "Delaunayhood" is stronger since not only does each triangle have an empty 
circumscribing ball but the ball can be chosen to be its diametric ball. Extending the notion 
of Gabriel graphs of a set of points in two dimensions, we call such a triangulation Gabriel. 

Theorem 2.2 For 6 = 2sin24e and e < any {£,6)-dense triangulation can be flipped to a 
Gabriel triangulation. 

Again, notice that the condition on uniformity becomes relaxed with increasing density. 
2.3 Dense triangulations 

It is natural to seek similar results for dense triangulations that are not necessarily uniform. 
It turns out that such triangulations can be flipped to almost Delaunay triangulations but not 
necessarily to Delaunay triangulations. To prove this result we will need some generalizations 
of the concept of stabbing as well as Delaunay triangulations. We denote a ball B = B{c,r) 
shrunk by a as B", that is, B" = B{c,r — a). With this definition, Df denotes the diametric 
ball of t shrunk by a. 

Definition 4 A surface triangulation T is a-Gabriel if for each triangle t £ T, the shrunk 
diametric ball Df contains no vertex of T inside. 



p 




B2 

Figure 1: (left) : pqr is stabbed by v and is locally stabbed by s; pq is a flippable edge, (right) 
: Bi and B2 are (/3)- and (— /?)-balls of pqr which is /3-stabbed by v and is locally /3-stabbed 
by s; is a /3-flippable edge. 
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Let nt denote the outward normal of a triangle t ^ T. For a triangle t ^ T and /3 G M, a 
/3-ball of t is a circumscribing ball of t whose center is at c + /Jn^ where c is the circumcenter of 
t. Observe that 0-ball of t is its diametric ball Dt- For any /? ^ 0, there are two balls of radius 
\/ p{tY + one is /?-ball and another is (— /3)-ball of t. See Figured) 

Definition 5 A triangle t (z T is /3-stabbed if a vertex of T stabs both (3- and {—(3)-balls of t. 
We say t is locally f3-stabbed if the stabbing vertex is one of the three neighbor vertices of t. 

Observe that if a triangle t is not /J-stabbed, the intersection of its /3- and (— /3)-balls cannot 
contain any vertex of T. This intersection contains the ball , a = p{t) + [3 — \J p{ty + 
which also cannot contain any vertex of T. Observe that a ^ (3. 

Observation 2.1 A surface triangulation T is [3-Gabriel if it does not have any (3-stabbed 
triangle. 

We prove the following results. 

Theorem 2.3 For e < 0.1, any e-dense triangulation of a surface with reach 7 contains a 
P-stabbed triangle only if it contains a locally (/? — 88e^^) -stabbed triangle. 

Choosing (3 = 88e^7 we conclude that there is a 88e^7-stabbed triangle only if there is a 
locally stabbed triangle. Therefore, if one gets rid of all locally stabbed triangles, there cannot 
be any 88e^7-stabbed triangles. In other words, the triangulation becomes 88e^7-Gabriel by 
Observation 12. 1[ 

Theorem 2.4 Any e-dense triangulation of a surface with reach 7 can be flipped to a 88e^7- 
Gabriel triangulation if £ < 0.1. 

2.4 Background results 

The following well known results on normal approximations will be useful in our analysis. 
Starting with work of Amenta and Bern [1], several versions of these results have been proved. 
We pick appropriate ones for our purpose. Let n^; denote the outward unit normal of S at a 
point X G S. 

Lemma 2.1 ([9l [3]) For any two points x and y inT, such that \\x — y\\ ^ for some e ^ 5, 
Zn^., ny ^ and Zn^, {y - x) ^ arccos(|). 

Following lemma is an oriented version of a result in [12]. Here we use the fact that T has 
a consistent orientation with E. Interestingly, this property does not automatically follow from 
triangles being small and vertex set being dense. 

Lemma 2.2 ( |12] ) Let pqr be a triangle in a e-dense triangulation of a surface. Assume thatp 
subtends a maximal angle in pqr. Then, fore < Zn^g^, % ^ arcsine+arcsin(-^(sin(2arcsine)). 

Combining Lemmas 12.11 and 12.21 one obtains the following corollary. 

Corollary 2.1 Let q be any vertex in a triangle pqr G T where T is a e-dense triangulation of 
a surface. Then, for e < 0.1, Znp^^, n^ ^ 7e. 

Define the dihedral angle between two adjacent triangles pqr and qrs as the angle between 
their oriented normals, that is, Znpqr,nqrs- An immediate result from Corollarv 12.11 is that the 
dihedral angle between adjacent triangles in a dense triangulation is small. 

Corollary 2.2 Let pqr and qrs be two triangles in a e-dense triangulation of a surface. Then, 
for e < 0.1, Zripqr, n^^s ^ 14e. 
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3 Flip algorithm 



Our flipping algorithm is very simple. Continue flipping as long as there is a flippable edge, 
that is, an edge incident to a stabbed triangle. 

MeshFlip(T) 

1. If there is a flippable edge e G T then flip e else output T; 

2. T := T' where T A T'; go to step 1. 

The first issue to be settled is the termination of Mesh FN p. It turns out that this simple flip 
algorithm terminates if T is a e-dense triangulation of a surface for e < 0.1. 

For convenience we introduce the notion of bisectors using power distance. The power 
distance pow(i?,x) of a point x € M'^ to a ball B = B{c,r) is ||c — x|p — r^. For two balls 
Bi and B2 in M^, the bisector C{Bi,B2) is the plane containing points with equal weighted 
distances to Bi and i?2- If Bi and B2 intersect, the bisector C{Bi,B2) is the plane containing 
the circle where the boundaries of Bi and B2 intersect. For two triangles pgr and pqs sharing an 
edge pq, we write Cpq = C{Dpqr, Dpgg). The following lemma establishes symmetry in stabbing. 

Lemma 3.1 Let pgr and pgs be two adjacent triangles where s stabs pgr. If Znpqr,npqs < ^, 
r stabs pgs. 

Proof. It can be shown that the bisector Cpq separates r and s if the planes of pgr and pgs make 
an angle larger than ^ or equivalently Znpqr,npqs < |. Let C^g be the half-space supported 
by Cpq and containing s. Clearly, Dpqs lies inside Dpqj- in Cpq as s is on the boundary of Dpgg- 
On the other half-space supported by Cpq which does not contain s, Dpqj- lies inside Dpgg- But 
this half-space contains r which is on the boundary of Dpgr- This means r is inside Dpgs- H 

If an edge incident to a stabbed triangle is flipped in a triangulation with dihedral angles 
less than ^, the circumradius of each new triangle becomes smaller than the circumradius of one 
of the two triangles destroyed by the flip. Actually, this is the key to prove that flip sequence 
to get rid of all flippable edges terminate. 

Lemma 3.2 Let T be a surface triangulation with dihedral angles smaller than ^. Let pg ^ T 
be an edge incident to a locally stabbed triangle pgr and pgs be the other triangle incident to 
pg. We have p{grs) ^ max{p(p(/r), /j(pq's)} and p{prs) ^ Tiiax.{p{pgr) , p{pgs)} . 

Proof. We prove the lemma for p{grs). The case for p{prs) can be proved similarly. Consider 
the bisector Cpq of Dpg,. and Dpgs, see Figure [21 Let Cpq be the half-space supported by Cpg 
containing s and C+ be the half-space supported by Crs containing p. 

By assumption the dihedral angle between pgr and pgs is at most ^. Then, Lemma 13.11 
applies to claim that r stabs pgs. 

Clearly the center of Dg^-s lies in the union Cg^ U Cgg. First, assume that contains the 
center of Dg^-s ■ Clearly Dg^s n is contained in Dpgr as s is contained in Dpg^ by the assump- 
tion that s stabs pgr. Therefore Dg^s is contained in Dpg^ in C+, which contains the center of 
Dgrs. This implies that Dg^s is smaller than Dpg^ establishing the claim. If Cgg contains the 
center of Dg^s the above argument can be repeated replacing Dpgr with Dpgg and s with r. H 

Since circumradii of triangles decrease by flipping flippable edges, triangles still can be ori- 
ented consistently with S and a homeomorphism using closest point map [2] can be established 
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Figure 2: Lemma l3.2t triangles after flipping have smaller circmnradii. 

between S and the new triangulation. In sum, the new triangulation satisfies the conditions 
for being e-dense. 

Corollary 3.1 IfT-^T' for a flippable edge e and T is e-dense, then T' is also e-dense. 
Lemma 3.3 IfT is e-dense for e < 0.1, MeshFlip(T) terminates. 

Proof. Let Ri, R2, .-, Rn be the decreasing sequence of the radii of the diametric balls of the 
triangles at any instant of the flip process. First of all, an edge flip preserves the number of 
triangles in the triangulation. An edge flip may change the entries in this sequence of radii, 
but not its length. We claim that after a flip the new radii sequence R'^, R2, R'^ decreases 
lexicographically, that is, there is a j such that Ri = R[ for all 1 ^ i ^ j and Rj+i > R'j+i- Let 
J + 1 be the first index where Rj+i 7^ ^'j+i- Since each flip maintains e-density (Corollarv l3.ip 
the dihedral angles between adjacent triangles remain at most 14e by Corollary 12.2 [ This angle 
is less than ^ for e < 0.1. One can apply Lemma 13.21 to each intermediate triangulation. By 
this lemma the maximum of the two radii before a flip decreases after the flip. This means 
the triangle corresponding to the radius Rj+i has been flipped and its place has been taken 
by a triangle whose circumradius is smaller than So the new radii sequence is smaller 

lexicographically. It follows that the same triangulation cannot appear twice during the flip 
sequence. As there are finitely many possible triangulations with a fixed number of vertices, 
the flip sequence must terminate. H 



4 Uniform dense triangulation 

We prove Theorems 12.11 and 12.21 now. First, we need some technical results (Lemmas 14. 2| 
and 14. 3p . We want to prove that if a vertex stabs the diametric ball of a triangle, it does not 
project orthogonally to a point inside that triangle. Next lemma is used to prove this fact. 

Lemma 4.1 Assume that a vertex v stabs a triangle pqr in a e-dense triangulation of a surface 
where e < 0.1. Let v be the point in pqr closest to v. The angle between the segment vv and 
the line of npgr is at least ^ — 26e. 

Proof. Let T be a e-dense triangulation of surface E with reach 7. Since v stabs Dpqj^^ we 
have \\p — v\\ ^ lej which implies that ||f — v\\ ^ 257. Walk from v towards v and let abc be 
the first triangle in T that we hit. Let y be the point in abc that we hit. (The triangle abc 
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could possibly be pgr.) We have ||f — y|| ^ \\v — v\\ ^ 2e^. By e-density assumption, we have 
||o — y\\ ^ 2ej. It follows that \\a — v\\ ^ \\a — y\\ + \\v — y\\ ^ 4e7. Then, Zn„,na ^ 8e by 
Lemma [27T1 and Zuabc'^-a ^ 7e by Corollarv 12.11 Therefore, Zny,iiabc ^ 8e + 7e ^ 15e. 

Let £ be an oriented line through v and v such that ^ enters the polyhedron bounded by T at 
y € a6c and then exits at v. Assume to the contrary that i makes an angle less than ^ — 26e with 
iipqr- Since — u|| ^ 2e7, Lemma [2T] and Corollarv 12.11 implv that Zni,,npgr ^ 4e + 7e ^ lie. 
Thus, i makes an angle less than ^ — 15e with n„. Since Zny,nabc < 15e, i must make an 
angle less than ^ with Uabc- Because i enters at y and then exits at v, Zn„, iiabc is greater than 
TT — (| — 15e) — f = 15e, contradicting the previous deduction that Zn^,nabc < 15e. H 



Lemma 4.2 Assume that a vertex v stabs Dpgr of a triangle pgr in a e-dense triangulation 
where e < 0.1. There exists an edge, say pq, such that r and v are separated by the plane Hpq 
that contains pq and is perpendicular to pqr. 

Proof. By Lemma |4.H vv makes a positive angle with the line of Upgr- It follows that v does 
not project orthogonally onto a point inside pqr. Hence, there exists an edge pg such that Hpq 
separates r and v. □ 

Next lemma leads to Theorem 12.11 This is where we require bounded aspect ratios of 
triangles which ultimately lead to the uniformity condition. The aspect ratio of a triangle t is 
the ratio of p{t) to its smallest edge length. 

Lemma 4.3 Assume that a vertex v stabs a triangle pqr in a e-dense triangulation T where 
each triangle has aspect ratio a < 2sin24E • V ^ < either pgr is locally stabbed or v stabs a 
triangle t such that pow{v, Df) < pow{v, Dpgr). 

Proof. By Lemma 14.21 there is a plane Hpq through the edge pg and perpendicular to pgr 
such that Hpq separates r and v. Let pqs be the other triangle incident to pg. If s lies inside 
Dpqr, pqr is locally stabbed and we are done. So assume that s does not lie inside Dpgr. By 
Corollary [22] Znpgr, ripgs ^ 14e, which is less than ^ for e < ^. Therefore, Hpg separates r 
and s too. It means that v and s lie on the same side of Hpg] see Figure [3j 

Let Cpg denote the bisector of Dpgr and Dpqs. Suppose that Cpq contains v and s on 
the same side. It follows that Dpqg contains Dpgr inside on this side as s lies outside Dpgr. 
Also V lies inside Dpgg since v lies inside Dpgr. It immediately implies that v stabs pqs and 
pow(v, Dpqs) < pow{v, Dpgr). Therefore, we can establish the lemma if we prove that Cpq 
contains v and s on the same side. This is exactly where we need bounded aspect ratios for 
triangles. 

Let s and v be the orthogonal projections of s and v respectively onto the line of pq. Consider 
the following facts. 

(i) The acute angle between ss and rip^j. is equal to | — Znpgr, n^^s, which is at least ^ — 14e 
by Corollary [221 

(ii) The angle between Hpg and Cpq cannot be larger than Zupgr, npgs which is at most 14e. 

(iii) We prove that /-Upqr^VV > Zllpgr, npgs — ZHpq, Cpq. 

The above three facts together imply that Cpg contains v and s on the same side as Hpq. 
Therefore, only thing remains to prove is fact (iii). 
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Figure 3: (left) : triangle pgr is stabbed by v. Both v and s lie on the same side of Hpg and Cpg. 
The case of v being in the thin wedge between Hpg and Cpg is eliminated if pgr has bounded 
aspect ratio, (middle) : the worst case for angle Zvpg. (right): the planes of Hpg and vpg make 
large angle ensuring v and s are on the same side of Cpg. 



First, observe that if v is the closest point of v in pq, we have by Lemma l4. II 

TT 



Znpgr,VV ^ 



26e ^ Ue ^ Zn 



pqr J Ilpgs ■ 



So, assume the contrary. In that case, the closest point of v in pq is one of p or q. Assume it 
to be p. Since v lies outside pg, the angle Zupg is obtuse. We claim that this angle cannot be 
arbitrarily close to vr. In fact, this angle cannot be more than the maximum obtuse angle pq 
makes with the tangent plane of Dpgr at p. Simple calculation (Figure middle)) shows that 
this angle is § + arccos \\p — g\\/2p{pqr) giving 

vr 1 

/-vpq ^ h arccos — 

2 2a 

where a is the aspect ratio of pqr. Since Dpgr contains v inside, \\v — p\\ ^ 2e'y. By Lemma |2. II 
ZupjVp ^ arccos e. Applying Corollary 12. II we get 

Zupgr, vp ^ Znp, vp — Zupgr, lip ^ arccos e — 7e. 

Let zpllf-u (Figure [3)[right)). Then, Zvv,vp = Zvpz = Zvpq — f ^ arccos One has 

1 



Zripqr, vv ^ Zupgr, vp — Zvv, vp = Ztipgr, vp — Zvpz ^ arccos e — 7e — arccos 



2a 



TT 



1 



^ lOe — arccos — for e < ^ 



2a 



72- 



We are now left to show that 



TT 



We — arccos ■ 



2a 



TT 



requirmg. 



24e > arccos ■ 



2a 



or, a < 



1 



2 sin 24e 

This is precisely the condition required by the lemma which can be achieved for e < 
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Proof. [Proof of Theorems 12.11 and 12. 2i ] The 'if part of Theorem 12.11 is obvious. To prove 
the 'only if part, let pgr be stabbed by v. With 6 = 2sin24e aspect ratios are at most 
1/(2 sin 24e). So, by Lemma 14.31 either pqr is locally stabbed or v stabs a triangle t where 
pow{v, Dt) < pow{v, Dpqr). In the latter case repeat the argument with t. We must reach a 
locally stabbed triangle since the power distance of v from the diametric balls cannot decrease 
indefinitely. For Theorem 12. 21 observe that maximum circumradius decreases after each flip and 
nearest neighbor distance cannot be decreased by flips. So, MeshFlip maintains (e, (5)-density 
after each flip which is the only thing remained to be proved. H 



5 Dense triangulations 

We establish Theorems 12.31 and 12.41 in this section. We drop the uniformity condition, i.e., we 
assume T is only e-dense for some e > 0. We will use the notation vrf to denote the plane 
parallel to t and passing through the point c + (3nt where c is the circumcenter of t. In other 
words, vrf is the diametric plane parallel to t in the /?-ball of t. 

Lemma 5.1 For e < 0.1, let T be a e- dense triangulation of a surface with reach 7. A vertex 
stabs a [3 -ball of a triangle t £ T only if there is a triangle t' € T with v as a neighbor vertex 
and V stabs the circumscribing ball oft' that has center in the plane o/vrf . 

Proof. Let i? be a /3-ball of t stabbed by w. Consider the edges of T lying in B and planes 
passing through these edges which are orthogonal to vrf . Let Pg denote such a plane passing 
through the edge e. Let t = pqr and pq be the edge so that Hpg separates r and w according 
to Lemma 14.21 The line segment rw must cross Ppg = Hpg and possibly others. Let pq = 
ei,e2,...,efc be the sequence of edges so that rw crosses Pe^^, Pe2, Pe^ in this order. 

Consider two triangles ti and fj+i incident to any edge Cj in the sequence ei, 62, e^. Let 
Bi and be the two balls circumscribing ti and tj+i respectively and having centers on the 
plane vrf. Observe that the bisector of Bi and -Bj+i is Pe^. If a vertex of tj+i lies inside Bi we 
have t' = ti satisfying the lemma. Otherwise, Bi is contained in -Bj+i on the side of Pg^ which 
contains w. So, pow(i?i, w) ^ pow(i?j+i, ty). Since this relation holds for any i G [1, k], we have 
either found the triangle t' satisfying the lemma or > pow(i?, w) = pow(Bi,w) ^ pow(i?fc, w). 
In the latter case pow{Bk,w) is negative and hence B^ contains w inside. The ball B^ circum- 
scribes tk and has center in vrf. It is stabbed by w where w is a neighbor vertex of t^ satisfying 
properties of t' required by lemma. H 

Proof. [Proof of Theorem 12. 3j Let t £ T he /3-stabbed. By definition, the /3-ball and (— /?)-ball 
of t are stabbed by a vertex w. Apply Lemma 15.11 to both of these balls. Observe that the 
planes (Pg) that we construct in the lemma remain same for both of these balls. It means 
that the segment rw in the proof crosses same set of planes. In other words, the triangle t' 
guaranteed by Lemma 15.11 remains same. Let B and B' be the two circumscribing balls of t' 
which have their centers in vrf and vr^^ respectively. If we prove that B and B' are larger than 
(/? - 88e27)-ban of t' , we will be done since then t' win be locally (/3 - 88e27)-stabbed. 

Let c and c' be the circumcenters of t and t' respectively. Since t' has an edge in the diametric 
ball Dt and all triangles have circumradius less than £7, the distance ||c — c'|| is at most 267. 
We have Znt,np ^ 7e and Znf/,ns ^ 7e where p and s are vertices of t and t' respectively. 
Also the distance between p and s cannot be more than 4^7 which gives Znp, ^ 8e. In all, 
Zrif, n^' ^ 22e when £ < 0.1. We want to estimate the distance of c' from the plane of t. In the 
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worst case this distance is 



||c — c'll sin Zrij, n^/ ^ Aej sin 22e ^ 886^7. 

It means if we choose /? > 88e^7, the center c' lies inside the slab made by offsetting vr^ by /3 
on both sides. The distance of c' from these planes is at least (/? — 886^7). Therefore, with 
(3 > 88e^7 we have B and B' larger than {(3 — 88e^7)-ball of t' proving the claim. □ 

Proof. [Proof of Theorem 12. 4j We apply algorithm Mesh Flip on the e-dense triangulation T of 
a surface whose reach is 7. According to Theorem 12.31 output triangulation cannot have any 
88e^7-stabbed triangle. By Observation 12.11 the output is 88e^7-Gabriel. □ 

Instead of flipping all locally stabbed triangles, one may flip more conservatively. If we 
go on flipping edges that are incident to /3-stabbed triangles, we get a triangulation which is 
(/3+88e^)-Gabriel. We flip less edges than MeshFlip does and hence obtain a worse triangulation 
in terms of approximation to Gabriel triangulation. 



6 Conclusions 

In this work we showed that a uniform dense surface triangulation can be flipped to a Delaunay 
one using simple Delaunay-like flips. If uniformity condition is dropped, we get almost Delaunay 
surface triangulation. 

This research ensues some open questions. Can the dense triangulations be flipped to exact 
Delaunay triangulation? It is unlikely that such triangulations can be flipped to exact Gabriel 
triangulation. It might very well be that they cannot be flipped to exact Delaunay triangula- 
tions. Our flip algorithm converts dense triangulations to almost Gabriel triangulations. Is it 
true that such triangulations are actually a weighted Delaunay triangulation of its vertex set 
weighted appropriately? Or, is it possible to assign weights to the vertices and carry out edge 
flips to convert a dense surface triangulation to a weighted Delaunay one? We plan to address 
these questions in future work. 
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